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Abstract: We study physical properties of strongly coupled massive quantum liquids 
from their spectral functions using the AdS/CFT correspondence. 

The generic model that we consider is dense, heavy fundamental matter coupled to SU(N C ) 
super Yang-Mills theory at finite temperature above the deconfinement phase transition 
but below the scale set by the baryon number density. In this setup, we study the current- 
current correlators of the baryon number density using new techniques that employ a 
scaling behavior in the dual geometry. 

Our results, the AC conductivity, the quasi-particle spectrum and the Drude-limit param- 
eters like the relaxation time are simple temperature-independent expressions that depend 
only on the mass-squared to density ratio and display a crossover between a baryon- and 
meson-dominated regime. We concentrated on the (2+l)-dimensional defect case, but in 
principle our results can also be generalized straightforwardly to other cases. 



Contents 



1. Introduction 



1 



2. Setup 



2 



3. Solving the Fluctuation Field 



5 



4. Physical results 

4.1 AC Conductivity 

4.2 Relaxation Time in the Drude limit 

4.3 Quasiparticle Spectrum 



8 

8 
10 
12 



4.3.1 Analytic Continuation 



12 



4.3.2 Results 



14 



5. Discussion 



17 



1. Introduction 

In the last decades, the AdS/CFT correspondence [1, 2, 3] has become a powerful tool 
to study various properties of the strong coupling limit of conformal field theories, with 
applications to QCD and also more recently to some aspects of condensed matter physics. 
In principle, one has to distinguish between top-down setups that are constructed within 
string theory and imply consistency and bottom-up setups in which the gravitational duals 
are constructed from a phenomenological point of view. In this paper, we use the former 
approach as we would like to explore what happens to a particular consistent theory. 

A common, and in particular top-down, approach how to introduce fundamental matter 
in AdS/CFT is the probe brane approach, where one considers a small number of "probe 
branes" in an AdS black hole background which are dual to the fundamental matter coupled 
to an adjoint gauge theory. For example one considers the well-known AdS§ x S 5 (black 
hole) solution (above the deconfinement phase transition) that is a solution to a stack 
N c 3> 1 D3-branes in the decoupling limit and is dual to a (thermal) SU(N C ) Af = 4 
supersymmetric Yang-Mills theory [2]. Then one inserts Nf <C N c intersecting Dp branes, 
giving Nf families of (charged) fields in the fundamental representation of the SU(N C ), 
living along the directions of the intersection [4, 5] - such as the above-mentioned D3-Dp 
intersections. 

In condensed matter applications, there has been particular interest in 2+1 dimensional 
systems that can be typically constructed using M2 branes [6], or as a defect in a 3 + 1 
dimensional background using D3-D5 [7, 8] and D3-D7 [7, 8, 9] intersections - and also 
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as bottom-up setups in various contexts such as superconductivity. As our world is 3 + 1 
dimensional, the defect setup may be more realistic, even though there are some problems 
with the consistency of the D3-D7 setup [7, 8]. In recent years, there has also arisen 
significant interest in quantum-liquid-like aspects that arise when the temperature is small 
compared to the density, for example in a 3 + 1 dimensional D3-D7 setup [10, 11], a 1 + 1 
dimensional D3-D3 configuration [12, 11] or in a 2 + 1 dimensional D3-D5 setup [11]. 
In these systems several physical aspects have been studied, such as the zero-sound mode 
[10, 13], a fermionic instability [14], the quantum Hall effect [15] and the response functions 
of fermionic matter at finite magnetic fields [16]. This low-temperature limit, in particular 
the construction of [11] will be in the focus of this paper 

One interesting property of these configurations is a scaling behavior that occurs at 
small temperatures in setups that are dual to a finite density of fundamental matter with fi- 
nite mass. Then one finds that solutions at different temperatures but fixed mass-squared to 
density ratios are equivalent through a simple scaling and display thermodynamic features 
reminiscent of quantum liquids [11]. In this paper, we will exploit this scaling to further 
explore the properties of this state of matter through their two-point functions. This limit 
in parameter space corresponds to large quark number density with finite temperature 
or finite quark number density with small temperature but still above the deconfinement 
temperature, and any choice of quark mass in each case. 

The paper is organized as follows: In section 2, we review the setup and the formalism 
developed in [11]. Then we will develop a method to solve the appropriate equations of 
motion in a temperature-independent parametrization in section 3 and obtain the physical 
results in section 4. There we first give an overview of the AC conductivity in section 
4.1, then obtain the small- frequency (Drude) limit in section 4.2 and the spectrum of 
quasiparticle excitations in section 4.3. Finally we conclude in section 5. 



2. Setup 

The supergravity background of a planar black hole in AdS§ is 

*°* = (" (l - fj <# + *%) + £ (j^p + *!§) , < 2 = (2.1) 

This corresponds to the decoupling limit of N c black D3-branes dual to M = 4 SU(N C ) 
super- Yang-Mills theory at finite temperature T, living along the flat directions of the AdS 
[3] . The temperature T is given by the Hawking temperature T = and the Yang-Mills 
coupling by gy M = 47rg s . Since the curvature L is given in terms of the string coupling g s 
and string length l s as L 4 = Air g s N c Z 4 , the 't Hooft coupling A = g Y yiN c can ^ e written as 

T 4 

A = -pf. Hence the "supergravity limit" L ^> l s in which the type IIB supergravity action 
and the solution (2.1) are valid corresponds to the strong coupling A S> 1. 
In practice, it is convenient to define dimensionless coordinates 

U:= 7' L2' x:=x Jj2> ( 2 - 2 ) 
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in which the metric g becomes 



ds 2 = ~ (-(1 - u^dt 2 + di 2 + + u 2 dn{\ . (2.3) 

In this setup all the fields transform in the adjoint representation of the SU(N C ). In 
QCD or condensed matter physics, however, one would like to consider also matter that 
is charged under this symmetry, i.e. that transforms in the fundamental representation. 
Hence, the use of this solution for such applications is very much limited. To introduce 
the fundamental matter one may create an intersection of "probe" Dp branes with the D3 
branes, such that at the string theory side there are fields at the massless level of field 
theory at the intersection. From the point of view of the probe branes, they correspond to 
endpoints of D3-Dp strings, and in the field theory dual they correspond to fundamental 
fields of a (defect) field theory. 

Here, we use the well-known D3-D5 defect setup (see e.g. [17, 18, 19]): 



12 3 

t x y z 

background : D3 x x x x 
probe : D5 x x x 



5 6 7 8 9 

9 (2-4) 



x x 



The dual field theory is now the SYM gauge theory coupled to Nf fundamental hyper- 
multiplets, which are confined to a (2+l)-dimensional defect. This construction is still 
supersymmetric, but the supersymmetry has been reduced from M = 4 to M = 2 by the 
introduction of the defect. In the limit Nf <C N c , the D5-branes may be treated as probes 
in the supergravity background, i.e. we may ignore their gravitational back-reaction. 

We assume translational invariance along the flat directions and rotational invariance 
on the sphere. Hence, the pullback on D5-brane gives us one scalar field corresponding to 
the position in the z direction, which was extensively studied in [7, 8], and another scalar 
which describes the size of the compact sphere and corresponds to turning on the mass 
of the fundamental matter, studied in [20, 8] and more extensively in the similar D3-D7 
system in [21, 22, 23, 24, 25]. Parametrizing the S 5 as rfQ§ = dO 2 + sin 2 6 dfl 2 . + cos 2 6 dQ?, 
and putting the branes on the first S 2 of the S 5 , the induced metric P[g] on the probe 
branes is given by 

ds 2 = J (-(1 - u^dt 2 + dS 2 + (-Lj + du 2 ) (2.5) 

+ l 2 (i - y 2 )dn 2 2 , 

where we defined the scalar as *S>(u) = sm(0(u)). 

We would like to turn on only the overall U (1) factor of the world-volume gauge field, 
so the the probe branes are governed by the DBI action 

S(G) = -T 5 N f [ V- det G , (2.6) 

J D5 
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where G = P[g] + 2irl 2 F. This U(l) gauge field is dual to the £7(1) current operator of the 
U(Nf) that gives rise to the "flavor symmetry". In particular, the flux of the electric field 

F = d u A t {u) du A dt , (2.7) 

corresponds in the field theory side to the baryon or quark number density (see e.g. [26]), 
where we explicitly keep the quark number density /0 quar k = iV c/ Obaryon 

P = ( Jt ) = -T^L = **Nf— KmduMu) . (2.8) 
5A t y r u->o 

From the action (2.6), it is straightforward to obtain the solution for the gauge field 

Mi « = pVXT ^_ , (2.9) 

^(1 - ^ 2 ){p 2 u A + (1 - * 2 ) 2 ) 



where p = „ . The equation of motion for \I/(u) is 




p2 u i + (1 _ ^2)2 



1 - * 2 )(1 - $ 2 + u 2 (l - u A ){d u ^) 2 
^(2(1- ^ 2 ) 3 - u 2 {d u t>) 2 (1 - u 4 ) (u A p 2 - (1 - ^ 2 ) 2 
u A J (1 - f 2 ) 3 (1 - ^ 2 + (d u V) 2 (1 - u 4 ) u 2 ) (u 4 p 2 + (1 - V 2 )' 



(2.10) 



+ 



On the horizon u = 1, this equation reduces to 

limd M ^ = -^- where * = lim (2.11) 

The asymptotic solution at u — > is 

$(n) = mu + cu 2 H , (2.12) 

where m and c are dimensionless parameters that are determined by the value of ^0 

on the horizon. Following arguments of the T-dual (3 + 1) dimensional D3-D7 setup 
[21, 22, 23, 24, 25], the quark mass M q and dual condensate C are given by 



M q = ^yTrh and C = T 2 N f N c c . (2.13) 

One can understand this identification of the mass from the separation between the D3 
and D5 branes in flat space, such that M q is the mass of a stretched D3-D5 string and the 
condensate is just the thermodynamic dual of the mass. 

At vanishing density, there is a critical temperature-mass ratio below which the probe 
branes do not extend down to the horizon [27, 21]. At finite densities, unless one turns on 
the scalar in the z direction considered in [8, 20], the branes always extend down to the 
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horizon even though a phase transition may still be observed at small densities [20] . In the 
limit considered in this paper, however, this phase transition is of no concern. 

To obtain the low-temperature scaling solutions, we make a transformation u = ^j\fp 
for (2.10) and expand it in the large-/5 (i.e. p <C T 2 ) limit. The leading order gives 



0, 



£ 2 y (1 - ^ 2 ) (1 - ^ 2 + (d^) 2 £ 



2^ 



(2 (1 - ^ 2 ) 3 _ (0 tf) 2^2 ^4 ^2)2> 

+ V j V JJ =0, (2.14) 



£ 4 (1 - J(l- * 2 ) (1 - * 2 + 2 £2) ^4 + (! _ ^2) 

which doesn't depend on p explicitly and there remains only implicit dependence from the 
scaling of £. In the large £ limit, we have a series expansion solution 

*«) = *»- *4=^ + -, (2.15) 

which serves as the boundary condition to solve (2.14) numerically. In the small £ limit, 
we obtain 

= ^£ + ce 2 + ... . (2.16) 
The parameters m and c are now related to the quark mass, condensate and density by 



m w /2iV c iV f , c 2C 

m=— = = M q \ — and c=- = — . (2.17) 

vp V /> A p p 



In general, equation (2.14) has no analytic solution and (2.15) implies that in practice 
one has to start integrating the equation from the large-£ region - that replaced the horizon 
boundary condition - to obtain the mass and the related condensate at the asymptotic 
boundary, rather than setting either of them in the beginning. 



3. Solving the Fluctuation Field 

To obtain the conductivity though linear response theory and also the spectrum of excita- 
tions, we need to know the two-point functions of current operators. In AdS/CFT, they 
are obtained from the variation of the on-shell action - since the baryon number current 
Ji is dual to the gauge field of the U(l) subgroup of the U(Nf), the variation with respect 
to Af 

/>2c: 

where dj(x - y) = -i0{x° - y°) ( [J^x), Jj(y)} > (3.2) 
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and Aifi is the boundary value of the gauge field at the asymptotic boundary u = 0. Hence 
we perform a pertubation F — > F + / of the gauge field and expand the action up to 
quadratic order, 

S = -T 5 N f [ y/-det(G + f) 
nS(G) + S®(G,f), (3.3) 

where 

S {2) (G,f) = ]r 5 N f [ V^detG^G^/^/^ - \g>°> f^G"? ' f aP ) , (3.4) 

4 JD5 1 

where G = P[g] + 1t:1 2 s F. The variation of eq. (3.1) then becomes 



Cij ~ ~ eo Mu) 



N f N c , . 

, £o = 2T^ 7 =^ (3.5) 



and we see that all we need to do is to solve the equations of motion for the gauge field 
obtained from eq. (3.4) in the background of sec. 2 which is all the usual well-known 
method. 

To obtain temperature-independent expression in the spirit of [11], we do this however 
in a smart manner that represents the temperature independence. To do so, we obtain a 
solution at large radius u <C 1 or £ <C \/p and fix it in an overlap in the region l/y/p <C 
u <C 1 or \fp 3> £ 3> 1, with an analytic solution in terms of the original radial coordinate 
u that is valid at relatively small radius 1/u <C \/p or £ 3> 1, and reflects the horizon 
boundary condition. 

We fix the gauge to A u = 0, and we work in the case of vanishing momentum where 
dxAn = and dyA^ = and the x and y component equations are the same. As discussed 
in [7, 8] in the context of the SX(2,Z) duality, the t-component is also related to the 
spatial components, such that we only need to obtain the solution for A%- After Fourier 
transformation, the x-component equation of motion is 



4, ((l-^ 2 ) 2 +nV)(l-^ 2 ) 

du I (1 - M) V (i-^ + (w^(i-« 4 )) ^ (<J ' u) 

+ (7^4) y (T^W) w u) = ' (3 - 6) 

where u) = -^L. In the near horizon region, (3.6) is reduced to 

ed e (ed e A s (u,u)) + j- A £ (oo,u) = (3.7) 
where e = 1 — u. The solution is 

M{u,u) = C + (l-u) + ^/ 4 + C-{l-uy iCj l\ (3.8) 
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where C± are integration constants. The in- falling boundary condition only picks up the 
second term. We set C_ to be 1 for convenience in the later discussion, which will not 
affect the conductivity and other physical results. 

First let us obtain the solution for 1/u <C yfp. In the limit p S> u~ 2 , (3.6) reduces to 

d u (u 2 (l - u A )d u A £ (uj, it)) + u 2 u) 2 A s (u;, u) + 0(p 2 u)- 2 = , (3.9) 

where we made use of the scaling in eq. (2.15). To solve this equation, we take the Ansatz 

Az(u,u) = CexpJ«" , (3.10) 

with some arbitrary number it*, the E.O.M. in this limit becomes 

u(-l + u 4 ) 2 ((u) 2 + (2 - 8u 4 + 6u s )((u) + uCo 2 + u(-l + u 4 ) 2 ('{u) = . (3.11) 

This equation can be solved order by order in the limit of large Co (i.e. 10 S> T) and we 
obtain the leading terms 

Hence the solution for A x (uj,u) is 

Ax(u}, u) = ao(u}) — exp(— itD-(ln[l — u] — ln[l + u] — 2 arctanktl)) , (3.13) 
u 4 

where ao(w) = exp(— iCo\ (ln[2] + f)) which is obtained by matching (3.8) and (3.13) in 
the near horizon limit. 

Next, we consider the E.O.M. in the outer region £ <C \fp- After transforming (3.6) 
into ^-coordinate, in the large p expansion the leading order gives the simplified E.O.M. 

(P(S)%M*> 0) + tfMfa + o(e 2 /p) 2 = o , (3.14) 

where 



a(0 



^ (1 - * 2 + (%#) 2 £ 2 ) (£ 4 + (1 - ^ 2 ) 

(e 4 + (i-* 2 ) 2 ) a-* 2 ) 



(3.15) 



0(0 = - , (3.16) 

and (D = tb/yfp. Note that the explicit p dependence has dropped out, which is the 
advantage of the ^-coordinate. In the large £ limit, this equation is further simplified to 

edt. (fdzM*, 0) + ^(i, = 0, (3.17) 

and the solution is 

A 5 (o;,0 = Ci-— + C 2 -— . (3.18) 
4 4 
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where C\ = and C2 = y / / 9e~ lWv/ ^ 3 ( ln ^ + 2 ) are fixed by matching this solution with (3.13) 
in the 1/y/p <n<l region. We can take (3.18) as the boundary condition to solve (3.14) 
numerically, replacing what would be usually the boundary condition on the horizon. 



Re A- x 



'/H 

- '/ \ 









0.0 , 0.2 0.4 0.6 0.8 1.0 



Figure 1: Comparison of the numerical solution for the full E.O.M. (3.6) (solid), the approximate 
solution (3.13) (dot-dashed) and the numerical solution of simplified E.O.M. (3.14) (dashed) with 
p = 70, Co = 50 (or Q = 5.98) and m = 34.1 (\I>o = 0.9, m = 4.08). The dotted vertical line indicates 
the location of 1/y/p. 



In figure 1, we see an example how the solution to the full equations, eq. (3.6), the 
"near horizon" solution (3.13) and the solution to (3.14) agree in the regions u 3> l/\/p 
and a < 1, respectively - and in the overlapping region 1/y/p <C u <C 1. 



4. Physical results 



In this section, we give an overview of some interesting properties that we can extract from 
the correlator eq. (3.1). For the transport properties, we obtain the conductivity from 
linear response theory, i.e. using the Kubo formula 

a%3 = -Cij . (4.1) 



4.1 AC Conductivity 

First let us look for an overview of the AC conductivity. Defining for convenience <7jj 
TrTaij, eqs. (4.1) and (3.5) give us 



a — —Eq 



^ = ^7=- (4.2) 



We should note that the conductivity doesn't explicitly depend on p, since it appears only 
in the pre- factor of Aj(£, Co) and gets canceled out in the conductivity. The AC conductivity 
for various in and Co is shown in the fig. 2. 
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5 10 15 20 25 30 

CO 

Figure 2: The real part of AC conductivity with various values of to. 

We see at large frequencies the "plasma" resonances that were discussed in [8] , and at 
small frequencies the absence of the Drude-peak. The latter is due to the fact that we took 
the limit w > T to obtain (3.13), and we will reconstruct it in the following section. We 
will also discuss the spectrum of "plasmons" that gives rise to the resonances in section 
4.3. 

An interesting observation is that the temperature-scaled AC conductivity goes to some 
finite value cfq (i.e. the conductivity scales as cto/(irT)) as the frequency Q approaches 
(but still above the Drude peak), which we plot in fig. 3 as a function of fa. 



0.500 - 




0.0 -^-> 0.001 — ' 1 ' 1 ' — ^ 

5 10 15 20 25 30 35 40 0.01 0.1 1 10 100 

m m 



Figure 3: The conductivity as a function of to. The dashed line indicates the asymptotic 
behavior Re(er) sa 8.94eo/w 2 - 

This constant, ^Re(c"o), approaches the constant 0.29 in the small fa region, and 
vanishes as Re(cio) ~ 8.94eo/m 2 when fa is large. Physically, it reflects an accumulated 
contribution from the Drude peak and all the plasmon modes, and the scaling implies 
that they contribute less as the quark mass becomes much larger than the scale set by 
the density - which is in line with the fact that the resonances become more stable, i.e. 
narrower with increasing mass as we will find out in section 4.3. 
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4.2 Relaxation Time in the Drude limit 

The Drude model describes usually the DC and small-frequency conductivity of weakly 
coupled systems, such as electrons in simple metals. In [8] it was found, however, that the 
Drude conductivity given by (see e.g. [28]) 

(4.3) 



1 — iujT 



also applies to strongly coupled systems and describes the conductivity very well in a 
consistent manner, up to a "conformal" contribution that is subleading at large densities. 
While the relaxation time was given only numerically, the DC conductivity was found to 
be 



-DC = ^V( 1 -^) 2 + P 2 -^- M) 

For convenience, we define the relaxation time f = T(Tiry/p) such that lot = Of. 

To obtain a closed- form expression for w < 1, we look for an approximate solution to 
(3.14) by considering only the first term. The solution for Ax(£) is then 

In the region <C £ <C 1/& region, /?(£) ~ £ 2 such that matching (4.5) with (3.18) 
in this region gives us C2 = 0. C\ cancels out in eq. (4.2), such that we obtain the 
conductivity 



a 1 

£ U) 



rHOdt') • (4.6) 



00 



This is purely imaginary and has a pole at oo = with the residue 

-1 



£0 



Res (-) = -( /" /rVX) . (4.7) 



Due to the Kramers-Kronig relations, the real part of conductivity becomes a <5-function 

in ui with a pre- factor 7rReso( — ) which we will compare to the Drude conductivity (4.3) 

£0 

by only considering the area under the peak, as our methods do not "resolve" its detailed 
frequency dependence. To do so, we integrate the real part of the conductivity eq. (4.3) 
over the frequency 

' " Re(a) du = f " du = (4.8) 

and identify this with the integral over the frequency in the low-frequency regime of our 
low-temperature-limit result, i.e. with the weight of the delta function. This is justified 
by the fact that our scaling and approximation oj 3> yfp has the effect of "shrinking" the 
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Drude-region to zero size in to. Using expression for the DC conductivity from [8] (4.4) and 
the expression for the residue, (4.7), we can obtain an expression for the relaxation time 



2p 



N c NfT 2 



+oo 



(4.9) 



At m = 0, this result can be solved analytically as /3 m =° \J\ + £ 4 such that 



2p 



+oo 



r(i/4) 2 



N c NfT 2 J TTTF 2 v^ N c N f T* 



(4.10) 




-0.04 -0.02 0.00 0.02 0.04 -0.04 -0.02 0.00 0.02 0.04 



Figure 4: The comparison of the conductivities obtained from a) the Drudc conductivity formula 
(4.3) (solid) b) our approximate small-temperature solution (4.6) (dashed) and c) the solution of 
the full equation (3.6), with p = 70 , rh = 4.08 and accordingly t = 196.72. 



In fig. 4, we compare the Drude conductivity with our numerically-obtained relaxation 
time to the result obtained from the full equations of motion and the result from our low- 
temperature limit for an illustrative large value p = 70 and finite value rh = 4.08. We 
see reasonably close agreement. Already for p = 100, there would be no difference visible 
between the result for the full conductivity and the one from the Drude conductivity. 



T 1 N c N f 




Figure 5: The relaxation time as a function of quark mass, m = ~?f y 2Jv ^ JVf 



region, r = tit 



2p 

N c N f 



approaches a constant, 
linearly with a slope of 0.27 when rh is large. 



r(i/4) 2 P 

2^N c N f T 2 ■ 



In the small m 
approximately 1.85, while it increases 
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In fig. 5, we show the relaxation time as a function of the quark mass. At small m, 
i.e. when the system is dominated by its density, we find r oc 5^2 W 1 From a weak- 

coupling "mean free path" point of view, this is somewhat counter intuitive, but it can 
be understood as the system is strongly coupled and the conductivity is due to collective 
excitations that become more stable as the density is increased. At large quark mass, we 
have t oc which is independent of the density as we are in the limit were the shortest 
length scale is the Compton wavelength of the quark mass (or more precisely of the energy 
scale of mesons [5] ^= that dominates the type of system at large quark mass [21, 24, 25, 8]). 
While this impression implies that the relaxation time increases with increasing mass and 
decreasing temperature, it is not of the form of the weak-coupling geometric scattering. 
Also the factor of the meson mass might be misleading, as they cannot carry the baryon 
number current. 

4.3 Quasiparticle Spectrum 

The spectrum of quasiparticles corresponds to the locations of the poles of the correlator in 
the complex frequency plane, and hence to quasinormal modes in the gravity side. Thus, it 
can be obtained by tuning the complex frequency ui such that the fluctuation fields vanish 
at the asymptotic boundary [29] . It is straightforward to search the exact locations of poles 
for large m by directly solving (3.14) for complex to numerically, but the numerics start 
to break down when we approach to the small m region, where we can only obtain few 
first order poles. Hence, to study for higher order excitations and small values of fh we 
developed a method based on analytic continuation. 

4.3.1 Analytic Continuation 

Our general strategy is to find a closed- form expression that is exact for large |u)| 3> 1. As 
this expression numerically still does not converge for sufficiently large imaginary parts of u), 
we then do a parametric fit for real frequencies and then perform an analytic continuation 
to obtain the location of the poles in the complex frequency plane. First, we take the 
Ansatz 

Mt) = Ceti«V#, (4.11) 

to re-cast (3.14) into 

«e) 2 + f|c(?) + CK) + ^ = o. (4.12) 

Then, we let £(£) = Co(£) + Ci(£)i where Co(£) is the leading term in u and satisfies the 
algebraic equation 

Co(0 2 + SfkoCO + = , (4-13) 



with the solutions 
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The subleading term Ci(£) should thus satisfy 

cue) + (2(0(0 + f|§) ci(o + c'oio + Ci(o 2 = . (4.i5) 

This is essentially a transformation of a linear second order equation into a non-linear first 
order one. 

The boundary condition at large £ can be obtained from the large-£ expansion of Ax(£) 
beyond the one of eq. (3.18), 



(4.16) 



which gives us 



Co(0 + Ci(0 = ^ HMO) « - £ + • • • ■ (4-17) 



By matching this with the large £ expansion of Co(£); 

& ^ w±to -r^-8^-i6^«---' (418) 

we know that we should choose the positive root for Co{0 and the resulting large-£ expan- 
sion of Ci(0 is given by 

which is taken as boundary condition for (4.15). It is further simplified as we can assume 
that Ci(0 is small compared to Co and /3' //3 such that (4.15) reduces to first order in our 
desired limit with the solution 

Ci(f) = < 2 ^ )+ H) (C C + jf* diA < 2 ^ )+i ^) C ' (o) , (4.20) 

where the integration constant is fixed by the boundary condition to = Ci(Ai) ~ 



* 4- 

rz + 



Given a solution for Ci, the solution for Aj(£) is then from the Ansatz (4.11) 



«<7 e j-|co(m'(i + / ci(0O, (4.21) 

J A 



where A is some large value of £. Our approximation holds for w> 1 when Ci(C') is 
small for some finite £. When £ approaches 0, the incoming and outgoing modes decouple 
such that (4.21) is still valid. This can be easily seen in eqs. (3.14) or (4.12) as the coupling 
is proportional to de («(£)/?(£)) which actually vanishes in both (large and small) limits of 
£. Thus, we could find the locations of poles by demanding J"_? Ci(£') d£' = —1. 
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In practice, however, we only solve (4.15) for real u, then we fit the real part of 
J+oo Ci(C') d£' into the real part of a generic complex model 



Re 



a e 



u c + f(u>) 



(4.22) 



where a, b and c are complex parameters, f(co) is a non-periodic function and asymptotically 
proportional to ui. Then, we can solve this model for Ci(£0 d£' = — 1 in the complex 
frequency plane to find the locations of poles in the complex frequency plane. One may 
obtain same result by fitting imaginary part accordingly and the result should also be 
independent of the particulars of the model. In particular, at large Co, the first term is the 
most generic one and as Co — > oo we are only interested in the parameter b to give us the 
tower of excitations. 

For validation of our procedure, we compare our results for each the first pole at various 
values of the mass in the range at which we can find the exact location to the exact result 
in fig. 6. Even though the first poles are expected to be the least reliable ones and most 
dependent on the model, they match with the exact locations very well. 
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Figure 6: The comparison between the exact locations of few first order poles (circle) and the 
approximate locations obtained from the analytic continuation method (square). The numbers are 
the corresponding values of in. 



4.3.2 Results 

Now that we have obtained the locations of the poles of the two-point function - exactly for 
large fh and with a good approximation also for small fh — we can use them to analyze the 
quasiparticle spectrum of our theory in the field theory side. To do so, we parametrize the 
poles corresponding to a fixed value of fh asymptotically (using the highest few excitations) 
as 

w n (m) = u!o(fh)n + P(fh) (4.23) 

for some complex parameters Coo(fh) and &(m) that represent the spacing of excitation 
levels (and their inverse lifetimes) and some "ground state energy" which will obviously 
not be exactly the value of the th excited state. The other way round, one can interpret 
the real part of l/ojo(m) as some induced length scale that depends on fh, i.e. on the ratio 
s/p/m. 
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That this parametrization is well-motivated can be seen in fig. 7, where we show the 
first few excitations for selected values of fh. 
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Figure 7: The locations of the first five orders of poles obtained from directly solving the E.O.M. 
for Aj(f) (3.14) (circles) and from the analytic continuation method (square) compared to our 
parametrization (4.23) shown as solid lines. The numbers in the graph are the corresponding values 
of fh. 



We see that parametrization is actually a very good fit, i.e. we have an almost linear 
tower of excitations, or an induced length scale that persists at high energies. We also find 
already now that both the excited states and the ground states get heavier and more stable 
with increasing quark mass. 

To look at the data from a different perspective, we also show the energy of the 
quasiparticle states, i.e. the real part of u> n , and their inverse lifetime to mass ratio as a 
function of the mass in fig. 8. 
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Figure 8: Left: The real parts of the locations of poles as functions of fh with linear fits shown 
as solid lines. Right: The double-logarithmic plot of — Im((D„)/Rc(w n ) as a function of fh. The 
solid line indicates the asymptotic behavior 1.12/ fh of the 5th order poles at large fh. The numbers 
indicate the order of the resonances and for the sake of clarity, the 3rd and 4th order poles which 
lie between the 5th and 2nd are not shown in the graph. 



We see that except for a small deviation at small mass (or large density), the energy 
of the excitations depends approximately linearly on the the quark mass. The inverse 
lifetime to mass ratio decreases inversely with increasing mass and, for the higher orders 



of the excitations, does not depend on the excitation number, such that the asymptotic 
behavior is approximately 1.72/m. 

This behavior can be best quantified in the parameters of eq. 4.23, that we show in 
fig. 9. 
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Figure 9: Top left: The "separation" of quasi-particles as a function of to, the solid line is the 
linear behavior obtained at large fh. The crosses correspond to values obtained with the fitting 
method explaind in the text. The solid and dashed lines show the limiting behavior in the large- 
and small mass regions. Top right: The "width to mass" ratio as a function of to, the solid line 
indicates the asymptotic behavior 2/to for large to. The plots on the bottom show the overall shift 
in the quasiparticles. To illustrate the level of reliability for the latter, we show the fits for the 
highest 2, 3 and 4 resonances. 

To explore the behavior or ujq at small masses that are out of reach of the method of 
sec. 4.3.1, we extract the data for J - ?^ Cl(£') fr° m the conductivity. From eq. (4.21), 
we get = Co(£) + Ci(0) an d since lim^o Co (0 = «I>, all the imaginary part of the 

conductivity originates from £i- As most of the contribution to the integral J £i comes 
from the region at small £ in which the two modes have decoupled already, we assume 
that a fit to lim£_>o (j (£) , i.e. to the imaginary part of the conductivity gives a good 
indication for the behavior of J £i and we perform a fit of the form HeAe 2 ™^^ to extract 
uiq. Interestingly, this method is reliable at small values of fh but not at large masses, and 
it is most reliable at fh = with loq ~ 2.20 + 2.20i. The energy gap between the excitations, 
Rea)o scales at large masses linearly with the mass as ujq ~ 12. Oy^r^r + ^-l^i. On the 
other hand, the overall shift D appears approximately constant, so it is not interesting for 
the rest of this discussion. Also, as we fit asymptotically, for the higher excitations, the 
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values are not reliable, as one can see from the dependence on the number of points take 
into the fit. The linear behavior with the factor of y/\ suggests that in the limit of large 
fh, energy levels of the quasiparticles are set by the meson mass scale that was found in 
[5], which suggests that the excitations, from fig. 7 we see even their ground states, either 
are now mesons, or at least that the there is scattering by mesons. For small mass, Oq 
approaches within the accuracy of the fittings the more precise value of u>q ~ 2.20 + 2.20i 
at fh = 0, such that the actual energy scale scales as ojq ~ 9.77y/^^y, so one of the 
open questions of [8] is actually resolved. Keeping in mind that p is actually the quark 
number density, this scaling suggests that the quark density arises mostly from baryons 
that again have plasmon excitations - density fluctuations that play a role in the context 
of condensed matter physics (see e.g. [30]) and are a pure quantum effect. As we already 
saw in fig. 8, the inverse lifetime to mass ratio starts with a finite value and decreases 
approximately as 1.95 2 where the prefactor depends slightly on how the date is fit and 
2 might be a limiting factor. The finiteness of this ratio just means that in all cases the 
lifetime decreases inversely with the excitation level. While in the plasmonic regime it is 
obviously independent from the quark mass, in the mesonic regime it increases proportional 
to fh which we can now interpret as the ratio of the meson mass to the plasmon mass, or 
the meson mass to the energy scale related to the geometric distance between the baryons. 
This also underlines the quantum nature of the decay of the mesonic excitations. We will 
discuss where this phenomenological appearance of baryons and mesons originates from in 
the conclusions. 



5. Discussion 

In this paper, we used the holographic technique to study various properties of fundamental 
matter with finite mass and baryon number density in the low-temperature regime, coupled 
to the usual N c S> 1,N = 4 SYM theory above the deconfinement phase transition. This 
state of matter might be similar in some aspects for example to what one finds in neutron 
stars. The fact that we used a top-down approach guarantees us that the field theory 
is consistent and well defined, giving more relevance to our findings. From the 2-point 
functions, we obtained on the one hand the AC conductivity over the whole frequency 
range in sec. 4.1, and on the other hand, we studied the underlying physical parameters 
- the relaxation time in the low- frequency Drude limit in sec. 4.2 and the quasiparticle 
spectrum and the nature of higher excitations in the large-frequency regime in 4.3. 

Using methods based on the approach in [11], we obtained results that are intrinsically 
temperature-independent except for the trivial dimensional temperature scaling. By a 
sensible choice of coordinates and taking the limit of a large quark number density to 
temperature-squared ratio in sections 2 and 3, we could eliminate the explicit density 
and temperature dependence from the equations of motion. Hence, the horizon, i.e. the 
temperature, appears only implicitly as zero-radius boundary condition and the density 
appears only implicitly through the density to mass-squared ratio. The relevant ratio is 
N P N so the relevant parameters are actually the baryon density and meson mass scale 
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rather than quark number density and quark mass. The role of the temperature as a 
dimensionful scaling parameter was replaced by the density. 

Looking first at the AC conductivity in sec. 4.1, we found instead of a drude peak a 
plateau that is increasingly suppressed with increasing mass, and then a set of resonances 
appearing with increasing mass before the conductivity reaches the conformal value [7, 
31, 32] at large frequencies. The plateau however appears not to be related to an energy 
gap due to the lack of an exponential suppression factor. Furthermore, from the Kramers- 
Kronig relations, we could extract and validate in sec. 4.3 the information of the Drude 
peak, that appears in our limit only as a delta function at zero frequency. We found 
an exact expression for the relaxation time that scales at small masses or large densities 
proportional to ^ \J jjjw^ so ^ actually increases with increasing density, which may be 
due to strong coupling as the fluctuations that carry the current are then more stable at 
large densities. At large quark mass or small density however the relaxation time scales as 
^^=, which is again not a classical mean free path type scaling. 

To obtain the underlying quasiparticle spectrum behind the resonances, we employed 
again a special scaling behavior to obtain a closed-form solution at large frequencies, and 
analytic continuation in sec. 4.3.1 to obtain the particle energies and lifetimes from the 
quasinormal modes in the gravity side. We saw in sec. 4.3.2 that we have an approximately 
linear tower of excitations with energy-independent width to mass ratio. At large mass, 
the energy scale of the excitations was set by the meson mass, whereas at small mass, it 
was set by N p Nf which could be interpreted as the energy due to the geometric length 
scale of the baryon number density. The inverse lifetime to mass ratio at large mass was 
approximately proportional to yj N p Nf and independent of the temperature. The latter 
implies that this lifetime is a pure quantum effect. At small mass however, it approaches 
a finite number. 

These scalings suggest that the quarks form at large density (small mass) baryons 
that have collective quantum excitations and at small density (large mass) the system is 
dominated by mesons. As discussed in [20, 11], we do not have phase transitions, but 
we are rather in the limits of an equilibrium between (not necessarily stable) baryons, 
mesons and free quarks. On the gravity side, this means that the embeddings of the probe 
branes converge towards the probe branes with Minkowski embeddings or probe branes 
of a baryon vertex (in a confining background), respectively. We illustrate this in fig. 
10, where we show the embeddings obtained in the small-temperature limit, eq. 2.14, for 
fh = 0.1 (solid black) and m = 10 (solid gray), scaled such that the axes correspond to 
r^se = Vp^I and nml = VIpl with mass oflfl = 3and ^ = 90 o and p = 0.09, 

respectively. For the small density, we compare it to the Minkowski embedding (black 
dotted) with p = and m = 3, that is the dual of the mesonic phase [5, 27]. The similarity 
in this regime is well-known [21, 20]. Stable baryons correspond to branes wrapping the 
KK wall in a confined background [33], forming a baryon vertex to which probe branes 
are attached and which source the electric field on the probe branes - in deconfining 
backgrounds like ours, a baryon vertex does not exist [34]. Now, the relevant scale is not 
the temperature (black hole radius ro) but the KK-scale, set by vkk, which represents the 
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Figure 10: The embeddings of the probe branes at to = 3 and m = 10 (p = 0.09, gray) and 
to = 0.1 (p — 900, black) compared to the Minkowski embedding of the mesonic phase (p = 0, 
dotted) and the baryon vertex embedding corresponding to m = 0.1 for a confinig background with 
Tkk of ro (dashed), ro/2 (dash-dot) and ro/4 (dash-dot-dot). The inset shows the baryon vertex 
and the KK-radius. 

deconnnement temperature. Following the instructions in [35] and [34], we show the baryon 
vertex embeddings corresponding to rh = 3 and p = 900 for confining backgrounds with 
tkk = ?~o (dashed), tkk = To/2 (dash-dot) and tkk = Tq/A (dash-dot-dot), motivated by 
the fact that we consider a deconfined background, i.e. we are above the deconfinement 
phase transition. We see that as tkk decreases, the embedding of the baryon probe brane 
approaches our embedding in the case of large densities - supporting the interpretation that 
the system is dominated by baryons in this limit. It is not quite clear, however, whether 
between these extreme limits there still free quarks, or whether the system consists entirely 
of baryons and mesons. 
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